Collective oscillations of a ID trapped Bose gas 
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Starting from the hydrodynamic equations of superflu- 
ids, we calculate the frequencies of the collective oscillations 
of a harmonically trapped Bose gas for various ID configu- 
rations. These include the mean field regime described by 
Gross-Pitaevskii theory and the beyond mean field regime at 
small densities described by Lieb-Liniger theory. The rele- 
vant combinations of the physical parameters governing the 
transition between the different regimes are discussed. 
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Recent experiments on trapped Bose gases at low tem- 
perature have pointed out the occurrence of characteris- 
tic ID features. These include deviations of the aspect ra- 
tio and of the release energy from the 3D behaviour 
as well as the appearence of thermal fluctuations of the 
phase, peculiar of ID configurations Interest in ID 
interacting Bose gases arises from the occurrence of quan- 
tum features which are not encountered in 2D and 3D. 
For example in ID the fluctuations of the phase of the 
order parameter rule out the occurrence of long range 
order even at zero temperature Such systems can- 
not be in general described using traditional mean field 
theories and require the development of a more advanced 
many-body approach. In the case of ID Bose gases in- 
teracting with repulsive zero-range forces, this has been 
implemented by Lieb and Liniger ||^ who studied both 
the equation of state and the spectrum of elementary ex- 
citations of a uniform gas. In the presence of harmonic 
trapping, ID Bose gases exhibit new interesting features. 
The corresponding equilibrium properties have been al- 
ready discussed in a recent series of theoretical papers 
(see and references therein). In the present work 

we investigate the consequences of harmonic trapping on 
the collective oscillations of an interacting ID Bose gas 
at zero temperature. We will consider various config- 
urations, ranging from the mean field regime [|| where 
the healing length is larger than the average interparticle 
distance, to the Tonks-Girardeau limit where the gas 
acquires Fermi like properties. We will show that the fre- 
quency of the lowest compression mode provides a useful 
indicator of the different regimes. 

We start our discussion from the hydrodynamic equa- 
tions of superfluids in ID 
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which decribe the dynamic behaviour of such systems 
at zero temperature. In these equations nxiz^t) is the 
ID density of the gas, v{z,t) is the velocity field, while 
Vextiz) is the external trapping potential which in the 
following will be assumed to be harmonic: Vext{z) = 
muj'^z'^/2. The hydrodynamic equations of superfluids 
have been already successfully employed to predict the 
collective frequencies of 3D trapped Bose-Einstein con- 
densates . A crucial ingredient of these equations is 
the local equilibrium (ie) chemical potential fige, which 
should be evaluated for a uniform ID gas {Vext = 0) at 
the density ni. The applicability of the above equations 
requires the validity of the local density approximation 
along the z-th direction. This is expected to be accu- 
rate for sufficiently large systems. Furthermore, Eqs.(|]j|) 
should be limited to the study of macroscopic phenom- 
ena where variations in space take place over distances 
larger than the average distance between particles. From 
Eq.(||) one can easily calculate the ground state profile 
through the equation 
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The collective oscillations are instead determined by writ- 
ing the density in the form ni(z, t) = ni(z)-|-e~*"*(Sni(z), 
with the function 5ni{z) obeying the linearized equation 
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which immediately follows from Eqs.(|^J|). In the case of 
uniform systems (Vext = 0) one has plane wave solutions 
Sni = eiip{iqz) with — c^q^, where q is the wave 
vector of the excitation and = ni{dfiie/dni)/m is the 
square of the sound velocity. It is worth reminding that 
the applicability of the hydrodynamic equations is not 
limited to the mean field scenario. Actually in ||^ it has 
been proven that in ID Bose gases the velocity of sound, 
derived from the macroscopic compressibility, coincides 
with the one derived from the microscopic calculation of 
the phonon excitation spectrum also in regimes far from 
mean field. 
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In the first part of the work we evaluate ni^ini) and 
the corresponding solutions for the collective oscillations 
in the framework of the mean field Gross-Pitaevskii the- 
ory. Even in this regime one can explore a rich variety 
of situations ranging from the Thomas-Fermi regime in 
the radial direction to the one of tight confinement where 
the motion in the radial direction is frozen. In the sec- 
ond part we extend the analysis to regimes beyond mean 
field, including the limit of the Tonks-Girardeau gas. 

Let us consider a uniform system of length L in the z-th 
direction and confined by a harmonic potential V{r±_) = 
in the radial direction. By writing the order 
parameter in the form = y^f{p±)/a±, where ni = 
N/L is the ID density, a± = y^h/mLUj_ is the oscillator 
length in the radial direction and p±_ = rj_/a± is the 
dimensionless radial coordinate, the 3D Gross-Pitaevskii 
equation yields the dimensionless equation 
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for the function / obeying the normalization condition 
27r / \f{pi_)\'^PA_dp±_ = 1. In Eq.(||), pge/hLU± is the chem- 
ical potential in units of the radial quantum oscillator 
energy. Eq. shows that the relevant dimensionless 
parameter of the problem is ani. It is worth considering 
two important limits. If ani ^ 1, one enters the radial 
Thomas- Fermi regime (hereafter called 3D cigar), where 
many configurations of the harmonic oscillator Hamilto- 
nian are excited in the radial direction and the equation 
of state takes the analytic form 
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Notice that in this limit the chemical potential is not 
linear in the density. This implies, in particular, that 
the sound velocity is related to the chemical potential by 
the law = /Ltfe/2m [|l^ rather then by the Bogoliubov 
relation = pie/m. A second important case is the 
perturbative regime where ani 1 (hereafter called ID 
mean field). In this case the solution of (||) approaches 
the Gaussian ground state of the radial harmonic oscil- 
lator and one finds the linear law 
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for the chemical potential. 

If the gas is harmonically trapped also along the z- 
th direction, the density ni exhibits a z-th dependence 
which is worth calculating as a function of the relevant 
parameters of the problem: the scattering length a, the 
number N of atoms and the trapping frequencies U!± and 
LUz- To this purpose one has to solve Eq.(^ by impos- 
ing the normalization condition J ni(z)dz = N. A useful 
quantity is the Thomas-Fermi radius Z defined by the 



value of z at which the equilibrium density ni(z) van- 
ishes. According to Eq.(|^), one has p, — pie{ani = 0) = 
(1/2)7710-1^^^. In terms of Z, Eq.(||) can be rewritten 
as flee{ani{z)) = {mojlZ'^ /2huj±){l — z'^/Z'^), where we 
have defined the dimensionless quantity /2£e(ani(z)) = 
[pie{ani(z)) — /i£e(arii — 0)] /huj±. This function is fixed 
by the solution of the Gross-Pitaevskii equation (^). Its 
inverse p,'^^ gives the value of art i as a a function of z, so 
that the normalization condition obeyed by the density 
can be written as 
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with t = z/Z. Eq. (g) explicitly points out the rele- 
vance of the dimensionless combination Naa± /a1 where 
= y^h/mujz is the oscillator length in the axial direc- 
tions. From Eq.(H) one can calculate, for a given choice 
of the parameters, the radius Z and hence the ID den- 
sity profile. In the 3D cigar limit Naa±/a1 ^ 1, one has 
Z = (a2/aj_)(15A^aa_L/a2)i/5 and 
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In the ID mean field limit Naa±/al <C 1, one instead 
finds M Z = (a2/a_L)(3iVaa±/a2)i/3 and 
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The density profiles are different in the two regimes, re- 
flecting the different behaviour of the equation of state. 
The conditions of applicability of the local density ap- 
proximation employed above are determined by requiring 
that Z 3> flz. In the ID mean field regime this implies 
the non trivial condition {az/a±){Naa±/a'^y^^ ^ 1. 

Let us now discuss the behaviour of the collective os- 
cillations. The hydrodynamic equation (^ has simple 
analytic solutions if the density derivative of the chemi- 
cal potential is a power law function: dpge/9ni cx nj~^. 
The Thomas- Fermi and the ID mean field regimes belong 
to this class of solutions with 7=1/2 and 7 = 1 respec- 
tively. Also the Tonks-Girardeau limit (see Eq.(|l3|) be- 
low) belongs to the same class with 7 = 2. Hence in these 
three relevant limits the dipersion relation for the collec- 
tive frequencies can be obtained analytically. By looking 
for solutions of the form nJ~^Sni{z) — z'^ + az'^^'^ . . ., 
where fc < 1 and only positive powers of z are included 
in the polynomial, one finds the result 
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The case fc = 1 corresponds to the center of mass motion 
whose frequency is given hy lo ~ ujz independent of the 
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value of 7. The most interesting k = 2 case (lowest com- 
pressional mode) is instead sensitive to the regime consid- 
ered. One finds uj^ = (5/2)u;^, lo^ = 3a;^ and uj^ — Auj^ 
for the 3D cigar, ID mean field and Tonks-Girardeau 
regimes respectively. The result = {5/2)uj'^ was first 
derived in [ pT[ by solving the hydrodynamic equations 
for a trapped 3D system in the limit of a highly elon- 
gated trap [llIz <JC uj±). This prediction has been con- 
firmed experimentally with high precision [ p^ . The re- 
sult uj"^ — Sw^ was derived in p^-p^, while the result 
— AujI for the Tonks-Girardeau gas is simply under- 
stood by recalling that, in this limit, there is an exact 
mapping with the ID ideal Fermi gas [p^ where the ex- 
citation spectrum, in the presence of harmonic confine- 
ment, is w = kujz- The same result has been recently de- 
rived inj 17| using the mean field equations of Kolomeisky 
et al. ilJl. 



of the harmonic confinement, so that the result for the 
collective frequency takes the compact form p9[| 
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FIG. 1. Transition between the ID mean field and the 3D 
cigar regimes: uj'^ /uJl as a function of the parameter Naa±/a^ 
(fuU line). Transition to the Tonks-Girardeau regime for 
a±/a = 100 (dashed line) and 10 (dotted line). 

In order to evaluate the collective frequencies in the 
intermediate regimes where the hydrodynamic equations 
are not analytically soluble, we have developed a sum 
rule approach [0, based on the evaluation of the ra- 
tio h'^uj'^ = mi/m^i between the energy weighted and 
inverse energy weighted sum rules. In the following 
we will limit the discussion to the lowest compres- 
sion mode which is naturally excited by the operator 
'^1=1 ^i- The energy weighted moment is given by 

where (z^) — J ni(z)z'^dz/N is the average square radius 
fixed by the ground state solution ni{z). The inverse en- 
ergy weighted moment m_i is related to the static polar- 
izability a by rri-i = (l/2)a. This is evaluated by adding 
the pertubation — ez^ to the Hamiltonian and calculating 
the corresponding changes 5{z'^) of the expectation value 
of the square radius: a = NS{z'^)/e. Adding the pertur- 
bation — ez^ is equivalent to changing the frequency lUz 



uj^ = -2 



d{z^)/dLjl 



(12) 



We have calculated the ratio w^/w^ as a function of 
the dimensionless parameter Naa±/a1, thereby explor- 
ing the transition between the 3D cigar and the ID mean 
field regimes. The results are reported in Fig. |^. For the 
experimental conditions of [|l]j2| where Naa±/a1 — 0.24 
and 0.08, we predict uj'^/ojI — 2.85 and 2.91 respectively, 
confirming that those experiments are actually touching 
the transition between the two mean field regimes. 

It is worth noticing that result ( |l^ ) for the collective 
frequency, being based on general sum rule arguments, 
applies also to regimes beyond mean field where the den- 
sity profile ni(z) and hence the value of (z^) cannot be 
evaluated starting from the Gross-Pitaevskii equation of 
state. Deviations from the mean field regime become im- 
portant when the healing length ^ — (87ran)~^/^ is com- 
parable to the average distance d between particles. In 
the presence of tight radial confinement one can use the 
relationship n = ni/i:a\ between the 3D density eval- 
uated at r_L = and the ID density ni = J n{r±)dr±. 
When a becomes smaller than d, one can write d = 1/ni . 
One then obtains the result ^/d — y/a^n^JSa which be- 
comes smaller and smaller as the ID density decreases, 
thereby suggesting the occurrence of important devia- 
tions from the mean field behaviour for very dilute ID 
samples. This should be contrasted with the 3D case 
where the mean field condition (^ > d) is better and bet- 
ter satisfied as the density decreases. The combination 
a'^ni/a can then be used as an indicator of the applicabil- 
ity of the mean field approach. When its value becomes 
of the order of 1 or smaller, one enters a new regime char- 
acterized by important quantum correlations. The corre- 
sponding many-body problem was investigated by Lieb 
and Liniger |^ who considered ID repulsive zero-range 
potentials of the form gioS{z). The interaction parame- 
ter giD can be also written as gm = fi.^ / (ma 1^) where 
aiD is the one-dimensional scattering length |2y]. By 
averaging the 3D interaction A'K^fi^a5{r)/m over the ra- 
dial density profile, one obtains the simple identification 
aiD ~ a\/a The comparison with the expression 

for ^/c? shows that the deviations from the mean field 
increase by decreasing amni. 

In the Lieb-Liniger scenario the energy per particle 
e(ni), when expressed in units of the energy /2mafjj, 
turns out to be a universal function of the dimensionless 
parameter aioni. Important limits are the high density 
limit amni 3> 1, where one finds e(ni) = h^ni/maio 
and hence the mean field result (^ for the chemical 
potential /i^e = d(nie{ni))/dni (a part from the con- 
stant term arising from the radial external force). Note 
that under the assumption a± ^ a, the two conditions 
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an-i <C 1 and amni » 1 required to realise the ID mean 
field regime can be simultaneously satisfied. The other 
important limit is the low density Tonks-Girardeau limit 
aioni ^ 1, where the chemical potential takes the value 
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Here the chemical potential no longer depends on the 
interaction coupling constant and reveals a typical Fermi- 
like behaviour 

In the presence of axial harmonic trapping the ground 
state density profile has been evaluated in Q using the 
local density approximation (^) . In this case the normal- 
ization condition J ni{z)dz = N takes the form 
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Similarly to Eq.(||), we have introduced the radius Z at 
which the density vanishes and the inverse of the function 
P-ieiniDam), where p,£e is now the chemical potential 
expressed in units of fi^ /2ma^jy. Eq.([l^) shows that the 
relevant combination of parameters is given by Na\j-,/o?^. 
This differs by the factor [a^/aY from the combination 
Naa\/o?^ characterizing the transition between the mean 
field regimes discussed in the first part of this work. 
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FIG. 2. Transition between the Tonks-Girardeau and the 
ID mean field regimes: uj^ /uj^ as a function of the parameter 
Nala/al. 

Analytic solutions are obtained in the two limits 
Nafjj/al > 1 and Nafjj/al <^ I. In the first case, 
one recovers the mean field result (|l0|). In the second 
one, we find the profile |p|,|l8| 
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with Z — \/2Naz. In this case the applicability of the 
local density approximation simply requires N ^ 1. 

By determining numerically the density profiles in the 
intermediate regimes we have calculated the frequency 



of the lowest compressional mode using the sum rule for- 
mula (^2|) . The results are reported in Fig. || as a function 
of Naljj/al. Fig. |^ shows the evolution of the collective 
frequency as a function of the parameter Naa±/a1 for 
two different choices of the ratio a±/a. The correspond- 
ing curves reveal the transition between 3D cigar, ID 
mean field and Tonks-Giradeau regimes. It is however 
worth pointing out that if the condition 3> a is not 
well satisfied, the ID mean field regime is not clearly 
identifiable (see for example the dotted line in Fig. |^). 
In this case the determination of the collective frequen- 
cies requires the simultaneous inclusion of 3D and beyond 
mean field effects. 
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